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Abstract 

Real scalar fields are known to fragment into spatially localized and long-lived solitons called 
oscillons or /-balls. We prove the adiabatic invariance of the oscillons//-balls for a potential 
that allows periodic motion even in the presence of non-negligible spatial gradient energy. We 
show that such potential is uniquely determined to be the quadratic one with a logarithmic 
correction, for which the oscillons//-balls are absolutely stable. For slightly different forms of 
the scalar potential dominated by the quadratic one, the oscillons//-balls are only quasi-stable, 
because the adiabatic charge is only approximately conserved. We check the conservation of 
the adiabatic charge of the /-balls in numerical simulation by slowly varying the coefficient of 
logarithmic corrections. This unambiguously shows that the longevity of oscillons//-balls is due 
to the adiabatic invariance. 
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I. INTRODUCTION 


Real scalar fields are known to fragment into spatially localized and long-lived solitons 
called oscillons [1, 2] or /-balls [3]. The oscillons//-balls are known to arise for various 
types of potentials such as a double-well potential [2, 4], the axion-like potential [5], or the 
inflaton potentials [6, 7]. The peculiarity of the oscillons//-balls is its longevity [8-10]. 
While the properties of the oscillons//-balls have been extensively studied from various 
aspects [8, 11-16], yet it is unclear what makes them so long-lived. This is in sharp 
contrast with the other types of solitons such as Q-balls [17]^ or topological defects [20, 21] 
whose stability is guaranteed by the conservation of the global U(l) charge Q or by their 
topological nature. 

It was suggested by Kasuya and the two of the present authors (MK and FT) in Ref. [3] 
that the longevity of oscillons is due to the adiabatic invariance. Such soliton that is long- 
lived due to the conservation of the adiabatic charge was named as the “/-ball”, so named 
because the adiabatic invariant is often represented by I [22], in much the same way as 
the Q-balls. The adiabatic current was shown to be conserved in a certain case where 
the spatial gradient energy is negligible. It was also argued that the /-ball conhguration 
is (quasi-)stable only if the scalar potential is dominated by the quadratic term. It is 
worth noting that, using the conservation of the adiabatic charge, the held conhguration 
inside the /-ball was estimated analytically, which showed a remarkable agreement with 
the numerical simulation. This observation provided a strong support for the conjecture 
that the longevity of the oscillons/I-balls is ensured by the conservation of the adiabatic 
charge. 

In this paper, as a further step in the direction of Ref. [3], we hrst give a rigorous 
proof that the adiabatic invariant in the classical mechanics can be naturally extended to 

^ In Ref. [18] it was claimed that the stability of /-balls can be understood by relating them to the 
corresponding Q-balls. To this end, they introduced U(l) breaking operators, which, however, spoil 
the stability of Q-balls as pointed out in Ref. [19]. Also, as we shall see later, the LR mass term 
potential (1) plays a special role in stabilizing the /-balls, which is hard to understand in terms of 
(/-balls. 
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a classical field theory and the adiabatic current is conserved for a scalar potential that 
allows periodic motion. In contrast to the previous work [3], this argument does not rely 
on the assumption that the spatial gradient energy is negligible. We then show that such 
scalar potential that allows periodic motion is uniquely determined to be the quadratic 
potential with a logarithmic correction like 


vw 


-m 
2 


, 2/2 


l-K\n 



( 1 ) 


where m is the mass parameter, and K is the coefficient of the logarithmic correction. 
For /-balls to be formed, K must be positive. Such a logarithmic correction often arises 
as a radiative correction in many examples, and it determines the strength of non-linear 
effects. In particular, the J-ball radius is determined by K (and m). If the scalar poten¬ 
tial is slightly deviated from the above form, the adiabatic charge is only approximately 
conserved, and so, the J-ball will split into smaller pieces in the end. We also perform 
numerical simulations to con&m that the adiabatic charge of the /-ball is indeed con¬ 
served. Specihcally we vary the value of K sufficiently slowly with time (adiabatically) 
and follow the evolution of the /-ball conhguration to see if their behavior agrees with 
the analytic solution based on the conservation of the adiabatic charge. Our analysis 
shows unambiguously that the stability (longevity) of the oscillons//-balls is due to the 
(approximate) conservation of the adiabatic charge. 

The organization of this paper is as follows. In sec. II, we hrst provide a proof of 
the adiabatic current conservation, and then determine the allowed form of the scalar 
potential. We derive the /-ball solution and study their properties both analytically and 
numerically in Sec. III. In Sec. IV we follow the evolution of the /-balls when a coefficient 
of the logarithmic potential K is varied with time and show that the adiabatic charge I 
is indeed conserved when the variation of K is adiabatic. The last section is devoted to 
discussion and conclusions. 
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II. ADIABATIC CURRENT CONSERVATION 


In this section we first give a proof of the adiabatic cnrrent conservation for a scalar 
potential that allows periodic motion. We then show that such a scalar potential is 
uniquely determined to be the quadratic term with a logarithmic correction. 

A. Proof 

Now let us show that the adiabatic invariant in the classical mechanics can be naturally 
extended to a scalar held theory. With the use of a constant of motion for a strictly peri¬ 
odic motion, we show that the adiabatic current is conserved while an external parameter 
is varied sufficiently slowly with time. Our argument and notation are based on Ref. [23]. 
We consider a scalar held theory with the following Lagrangian: 

C = U,4>d^ct>-V{ct>,a{t/T)), ( 2 ) 

where a{t/T) is an external parameter that varies sufficiently slowly compared to the 
typical time scale of the scalar held dynamics. T dehnes the time scale over which the 
external parameter a{t/T) changes from a* = a(0) to aj = a(l), and it will be set to be 
inhnity in the end. The Hamiltonian density is given by 

= + + V{(j),a{t/T)), (3) 

and the Euler-Lagrange equation is 

^{x) - + R'(0(x), a{t/T)) = 0, (4) 

where the prime represents a partial derivative with respect to (p{x), and the overdot 
means the time derivative. Using the equation of motion, one can write down a (non- 
)conservation law of the energy: 

_ dVda{t/T) _ mda{t/T) 

~~d^ dt ' 

where 

f = n, f = -4){x)di(j){x), ( 6 ) 
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where di = djdx^. The energy is not conserved because of the external parameter a{t/T). 
As is clear from the derivation, spatial components of the current arises from the gradient 
term, which is the crucial difference from the case of the single degree of freedom in 
classical mechanics. 

For later use, let us rewrite the above equation as 




- di 


d(j){x) 

dt 



dH da{t/T) 
da dt 


( 7 ) 


One can dehne another energy density "H which differs from "H by a total spatial derivative 
as 

n = n- ^di{ 4 >di(j)). ( 8 ) 

For a vanishing surface term, the spatial integrals of H and "H are equal: 


/ 


d^xH 


/ 


d^xH. 


(9) 


Using H, one can rewrite Eq. (7) as 

dU da{t/T) _ did da{t/T) 
da dt da dt 

This equation will be important in the following argument. 

We limit ourselves to the case in which the scalar dynamics is approximately periodic. 

In particular, we assume that, if a{t/T) is hxed to be constant, i.e. a{t/T) = a^r/T), the 

scalar dynamics is exactly periodic and the scalar held has a solution in a separable form. 




- {4){x)di(j){x) - ((){x)di(}){x] 


0(x) = <F(T, a(r/T))/(t, a{T/T)), (11) 

where f{t,a{T/T)) is a periodic function: 

/(t,a(r/T)) = f{t + l/ur,a{T/T)), (12) 

where Ur is the frequency of the scalar dynamics for a = a{T/T), and the maximum value 
of fit^aij/T)) is normalized to be unity. We emphasize here that such periodic motion 
is not guaranteed at all for a generic form of the scalar potential, and the scalar potential 
must be close to the quadratic one, as we shall see later in this section. Here we do not 
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specify the form of the potential in order to include a case in which the scalar dynamics 
can be approximated by the above separable form over a sufficiently long time scale of 
interest. Most importantly, H is a constant of motion for the separable solution (11) with 
a constant a = a^r/T), because <j){x)di(j){x) = (f){x)di4>{x) holds in this case (see Eq. (10)). 

As the external parameter a{t) depends on time, the scalar dynamics is not strictly 
periodic. In particular, a signihcant amount of the energy can be transferred to other 
spatial points by scalar waves, in contrast to the case of one dynamical degree of freedom 
in classical mechanics. For a constant a(t), the trajectory of (0(a:), 7r(x)) in the phase 
space is a closed one so that the modihed energy density Ti at each spatial point is a 
constant of motion. Here 7r{x) = (j){x) is the conjugate momentum. 

Following the approach of Ref. [23], we consider a hypothetical system for one period 
of the motion from t = T to t = T + l/h'r while the external parameter a{t/T) is hxed to 
be the value at f = r, i.e., a(t/T) = aij/T). In such a hypothetical system, the trajectory 
is a closed one, and this is indeed possible for a certain class of the scalar potential. We 
denote the trajectory of (0(x), 7r(a;)) in such a hypothetical system by (0T-(a^), 7rT-(x)). As 
mentioned before, for the separable form (11), one can End a constant of motion at each 
spatial point, 

H(7r^(x), 0^(x), a(r/T)) = p(r, x). (13) 

Solving this equation for n^^x), we can express vrT-(a;) as 


Tir = 7r^(0^(a:),p(r,f),a(r/T)). 


(14) 


So, tTt- can be regarded as a function of 0r(a;), p{t,x), and a^r/T).^ For later use, let us 


differentiate Eq. (13) with respect to p and a: 

/ d'KT-\ ( dH \ 

\dp)4>r,a 


(15) 

(16) 


^ The modified energy density depends on the spatial derivatives of (j){x). One may explicitly show such 
spatial derivatives, which, however, does not affect the following arguments. This is because what we 
need is the partial derivative of R or H with respect to tt and a. 
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where we used the following relations, 


m 

d^T 

m 

da 


(pT 


dn 

dn^ 


(dn 

\ da 




( 17 ) 

( 18 ) 


In analogy with the argument in classical mechanics, let us estimate the area in the 
phase space surrounded by the trajectory at each spatial point, which is given by 


rvr 

= 2 / 7r^(0^(x),p(r,f),a(r/T))(i(/)^, 


T+I/Vt 


(Prit, x)Ydt 


(19) 

( 20 ) 


where and are the two roots of vr^ = 0 and we assume 0^^^ < (0^^. Note that x 
represents {t,x) not (r, if) here. Let us hrst differentiate J^{t,x) with respect to r: 


dJ0x) ^ -7 / 

= 2 [7r^(0^(x),p(r,)r),a(r/T))]Li) 


dr 


1,(2) 


1 ,( 1 ) 


dl^r 

dp 


dp{T, x) ^ f dn. 


dr 


da 


da{T/T)) 
dr 


d0r, (21) 


'4>t' \ \ ~r / (p^^a ~ ■ \ / <pT,p 

where the first term vanishes as nj. vanishes at the end points, and note that 0,- in the 
integrand is an integration variable. Using Eqs. (15) and (16), we obtain 


dJ^{T,x) 0^ (dp{T,x) f da{T/T)) 

dr \ dr \da) ^ j, dr 


4>t'\ V / 7rT,0T j (f>r,a 


dn 


j 


-1 


r+l/Ur 


d pjr, x) 
dr 


dn 

da 


da^r/T)) 


TTt,? 


dr 


dt, 


where we have used 


dn \ 

dn.,) 




00r(x) 

dt 


#r, (22) 

(23) 

(24) 


^ We limit ourselves to a simple case in which the trajectory in the phase space has only two roots of n, = 
0 (like an ellipse). The extension to a more complicated (but periodic) trajectory is straightforward. 
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in the second equality. So far, there is no difference from the argument in classical 
mechanics in Ref [23] except for the extra label, x. We replace the time variable t in 
Eq. (10) with r, and substitute it into the above equation. 


dJ^{T, x) 

Wr 


r+l/i/r 


da{T/T) ) 

dr 

( (9'R(7r^(t,x),0^(t,x),a(r/T))^ 

da{T/T)) 

V daijjT) y 

dr 


dt 


+ a I 2 


r-r+l/,., £) d't>(T,x) 

-di(t){T,x) - (t)(T,x)di- 


dr 


dr 


dt 


( 25 ) 


where it should be noted that the third integral over t is trivial and the integrand is 
independent of t. Let us now dehne the spatial component of the adiabatic current as 


J^{t,x) = 






dr 


dr 


(26) 


and then, Eq. (25) can be rewritten as 


5^J^(r,f) = 


fT+l/l/r 


/5'H(7r(r,x),0(r,x),a(r/T))\ da{T/T)) 

\ daij/T) J dr 

d'H{7iT-(t, x), (j)r(t, x), a(r/T)) \ daij/T)) 


dt. 


( 27 ) 


da{T/T) J dr 

In order to show the conservation of the adiabatic charge, let us integrate (27) over one 
period from r = r* to r = r* + l/z/^-p 


>'Ti + l/Ur, 


rr+l/iXr 


^ fda{T/T))\ 

dTd^j!^{T,x) = dr ( - — - j 

g'R(7r(r,x),0(r,f),a(r/T)) \ _ / g'R(7r^(t, x), 0^(t, x), air/T)) 


da{T/T) 


da{T/T) 


(28) 


Now let us see that the above quantity approaches zero faster than C>(l/T) as T —)■ oo, 
which is necessary to show the adiabatic charge conservation (33). One can see that 
the RHS of Eq. (28) contains a factor, da(r/T)/(ir, which is proportional to 1/T. In 
addition, as we shall see below, the hrst and second terms contain an additional factor 
which oscillates fast about zero; the hrst integrand in the RHS is independent of t, and 
it contains functions 7r(r, x), and 0(r, x), which oscillate fast as r varies. In general it 





























oscillates fast about some finite value. The second integrand exactly subtracts the hnite 
value, as it is obtained by averaging the hrst term over one period. (Note that, in the 
limit of T —)■ oo, the different between 0 (tt) and 0,- (tt,-) becomes negligible.) Thus, when 
integrated over the period, the sum of the first and second terms approaches zero faster 
than C>(l/T). 

To summarize, we have proved that the adiabatic current is conserved, 

^ = 0 (29) 


with 

jo = — ^, r = -- Ud,<p - , (30) 

if the scalar held dynamics is periodic at each spatial point and if the external param¬ 
eter varies sufficiently slowly. Here u = 27iu is the angular frequency, and the overline 
represents the average over one period of the motion, i.e., 

pt+l/u 

X{t) = u X{t')dt'. (31) 

Note that the spatial components of the adiabatic current are induced by the weak de¬ 
viation from the separable form. This implies that the adiabatic charge is transferred to 
other spatial points gradually as the external parameter varies adiabatically, which allows 
deformation of the oscillons/I-balls as we shall see later. 

We dehne the adiabatic charge / as 


I = 



(32) 


where D denotes the spatial dimension and the pre-factor l/27r is just a convention. For 
a spatially localized conhguration, the adiabatic charge is conserved. 


/ = const. 


(33) 


as long as the external parameter changes sufficiently slowly with time. 
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B. Form of the scalar potential that allows periodic motion 


So far we have assumed the existence of a scalar potential V (0) that allows periodic 
motion for which the solution is given in a separable form, 


tp = $(£)/((), 


(34) 


where the periodic function f{t) is normalized so that its maximum value is equal to 
unity. Now we determine the form of such potential. Substituting the above separable 
solution into the equation of motion, we obtain 

f r($/) 

/ $ $/ • ^ ^ 
This equation implies that the derivative of the potential in the RHS should take a form 
of 




= kl($) + R(/), 


(36) 


where ^(*h) and B{f) are some functions of $ and /, respectively. On the other hand, as 
the potential V (0) is a function of 0, the derivative of the potential is given by 

r($/) 




= W), 


(37) 


where 0(0) is a function of 0. Combining the relations (36) and (37), we obtain the 
algebraic equation for C: 

C{^f) = kl($) + B{f) = 0($) + C{f) - kl(l) - B{1) 

= 0($) + 0(/)-0(l). 

Eq. (38) is satished if and only if C{(j)) = aln(0/6), where a and b are constants.^ Then 
we obtain 


(38) 


1 /( 0 ) = 


iCln 


2M2 


(40) 


This can be seen by noting that one can derive the following differential equation for C(0), 

dC(0) _ C(0 + A0)-C(0) _ C(l + A0/0)-C(l) l 

d0 A0 A<f>^o A0/0 0 0 


(39) 
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where m and K are constants. Therefore, the scalar potential must be the quadratic 
potential with a logarithmic correction, and we call it as the logarithmically running (LR) 
mass term potential in the following. Note that there are in fact only two independent 
parameters, as M can be absorbed by rescaling m and K. The above argument does not 
hx the magnitude and sign of the parameters. As we shall see in the next section, the 
J-ball solution exists if m? > 0 and 0 < A" 1. 


III. /-BALL SOLUTION 

In this section, we derive the /-ball conhguration as the lowest energy state for a given 
value of the adiabatic charge, using the LR mass term potential (40) with m > 0 and 
0 < /L -C 1. We will show that the /-ball conhguration is given by a Gaussian distribution, 
and we numerically conhrm that the scalar dynamics is periodic and p is a constant of 
motion for the /-ball solution. Note that the proof of the adiabatic current conservation 
and the form of the scalar potential are valid for any number of spatial dimensions D, 
and we consider the case of i/ = 1, 2, 3 in numerical simulations. 


A. Gaussian field configuration 


We would like to hnd a scalar held conhguration that minimizes the energy for a given 
adiabatic charge I in the same way as in the case of Q-balls. Using the method of the 
Lagrange multipliers, this problem is formulated as hnding a spatially localized solution 
which minimizes the following 


Ex = I d^x p{x) -|- A I / — / d^x — 


= XI + I d^x 


X\ 1- 


2 + 


(41) 


® In contrast to Q-balls, the conservation of the adiabatic charge I is ensured only for adiabatic processes. 
Therefore, our argument does not preclude the existence of configurations with a lower energy, which 
cannot be reached via the adiabatic process starting from our Gaussian solution. 
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where we have used p = p in the second equality as p is a constant of motion. 

For the separable form (34), one can perform time averaging of f{t). If there were not 
for the logarithmic correction, the periodic motion is simply given by a homogeneous scalar 
held oscillating in a quadratic potential. In this case, the periodic function is given by 
f{t) = cos(mf), and the time average of the oscillating functions is trivial: f{ty = 1/2, 
f{ty = w? 12, and f{ty\nf(ty = 1/2 — In 2. With the logarithmic correction, those 
results are modihed by a factor of 1 + 0{K), and we write them as 

JW = c, (42) 

/(f)2 = da;^ (43) 

/(t)2 1n/(t)2 = i, (44) 


where c, d and £ are constants of order unity. Then Ex is given by 

Ex = XI+ f d^x 


2 2 V oj J 


+ \ l--K\ -K\n 


(—) 

I 2M2 i 


(45) 


The bounce equation is obtained by taking a functional derivative of Ex with respect to 






2M2 


= 0 , 


where we have dehned as 


1 - K - -K + 1 - 2- 

C ) C \ OJ 


,d 


A 


(46) 


(47) 


We assume that the bounce solution is spherically symmetric, = 4)(r), where r is the 
radial coordinate. Then the Laplacian can be written as 


d^ 


D — 1 d 


dr^ r dr 




(48) 


Let us adopt the Gaussian ansatz [3] 


<F(r) = $cexp(-r^/i?^). 


(49) 
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where <hc is the amplitude of the /-ball at the center and R is the radius. Substituting 
(49) into the bounce equation (46), we obtain the relation as 

o 4 


/?2 


r 1 K j 


[2D , 9 

2 , 

— 



= 0 . 


(50) 


This relation (50) should be satished for an arbitrary value of r, thus the radius R and 
the Lagrange multiplier A are determined as 




(51) 


and 


UJ 


^=2 


c m 


1 + -,^{1 + {D-1)K--K-K\n 




2M2 


(52) 


From eq. (51), we can see that the radius of the J-ball is determined by the coefficient K 
and mass m. As mentioned before, the choice of M is arbitrary as it can be absorbed by 
rescaling m and K. If we set M = <Fc/\/2, the Lagrange multiplier is given by 




c m 




^u(,l + 0{K)), 


(53) 


where we have used a; ~ m and c ~ d in the second equality. 

Let us evaluate the adiabatic charge I for the J-ball prohle derived above, 

/ = / d^x—(f)^ = —/2 f dx'^d’^ exp(—2r^//?^) 


1 

m\l K 


vr 


u 


(54) 


where we have used (51). When the parameters are varied adiabatically, the adiabatic 
charge I is expected to be conserved. We shall see that this is the case in numerical 
simulations. 

For the Gaussian prohle, the modihed energy density p is given by 


P = ^0^ - 


= 


1 5. 


1 


-f + -m^f (1 + KD - K\n 


V2M2 J 


f- 


V 2M2 




(55) 
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Substituting the Gaussian solution (49) with (51) into the equation of motion (35), we 
hud that the periodic f{t) satishes 


/(t) +m^ 


1 + K{D -1) - K\n 


(- 


$2 


c 1-2 


V2M2 

This can be integrated to obtain the following relation, 


m 


fit) = 0. 


= m^{l - f) 


1 + KD-K\n 




2M2 


+ Krn^f^ In / 


(56) 


(57) 


where we have used the normalization, /(O) = 1 when /(O) = 0. Using (57) one can 
rewrite p as 

p{r) = pce~‘^''^tR^ ("53) 


with 


Pc 


2 


1 + DK -K\n 



( 59 ) 


For comparison with numerical simulations, we dehne the /-ball radius R 1/2 where the 
modihed energy density p is equal to pcl‘2^'- 


fhf2 1 , , 

We also dehne the effective amplitude of the scalar held, <hc in terms of the modihed 
energy density. 


<Fc = \/2pJrn? 

= + ( 61 ) 

Note that <hc is roughly equal to the actual oscillation amplitude <hc up to a correction of 
order K. 


B. Numerical simulations 

Here we numerically conhrm that the J-ball solution obtained above is indeed a solution 
of the equation of motion. In particular we will see that the modihed energy density p is 
a constant of motion. 
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The LR mass term potential (40) contains a logarithmic function of 0, and so we have 
inserted a small parameter e into the potential and its derivative as 




l-RTln e + 


2M2 


50 


= m 


1-K- 


-Kln e + 


2M2 


(62) 


e + ,p^/i2AP) 2AP 

for numerical stability. We have set e = 10“^° in our numerical simulations, and we have 
checked that our results are insensitive to the values of e as long as it is much smaller than 
unity. This regularization is adopted in the numerical simulations here and in Sec. IV. 

We have performed lattice simulations for the cases oi D = 1,2 and 3. As the initial 
condition we adopt the Gaussian prohle (49) with = 2M and K = 10“^, and followed 
its evolution from t = D to t = 10^/m. The box size L and the number of grids N for 
D = 1, 2, and 3 are 


' N = 2048, L = 100/m, for D = l 

< N = 256^ L = 100/m, for D = 2 (63) 

^V = 128^ L = 50/m, for D = 3 

for which the spatial resolution is Ax = 4.8 x 10“^/m, 0.39/m and 0.39/m, respectively. 
We set the time step as At = 10“^/m. 

We show the results for H = 1 in Fig. 1. In the top left panel, the spatial distributions 
of the modihed energy density pc at different times are shown. All the lines are overlapped, 
implying that the Gaussian ansatz is valid and pc stays a constant in time. From the other 
panels, we can see that all of pc, R 1/2 and <I>c remain constant in time and their values 
are in a perfect agreement with the analytic results (59), (60), and (61), respectively. We 
have also conhrmed that the time evolution and the properties of the J-ball conhguration 
in numerical simulations are in a very good agreement with the analytic results for the 
case oi D = 2 and D = 3. Therefore we conclude that the adiabatic charge I is indeed 
conserved in the numerical simulations. 
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FIG. 1: Numerical results of the I-ball for D = \. We have set K = 10“^ and <hc = 2M. The 
top left panel shows snapshots of the spatial distribution of p at t = 0,10,10^, 10^[l/m], and all 
the lines are overlapped, implying that p is a constant of motion. The top right, bottom left 
and bottom right panels show the time evolutions of pc, R 1/2 *^*0 in a- very good agreement 

with the analytic results (59), (60), and (61), respectively. The box size is L = 100/m and the 
grid number is N = 2048. 


IV. ADIABATIC DEFORMATION OF /-BALLS 

In the previous section we have derived the J-ball solution so that it minimizes the 
energy for a given adiabatic charge I in much the same way as Q-balls. We have also 
numerically conhrmed that the obtained /-ball solution indeed satishes the equation of 
motion and the modihed energy density remains a constant of motion, which plays a 
crucial role in the proof of the adiabatic current conservation. In this section, in order to 
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further support the conjecture that the stability (longevity) of the oscillons//-balls is due 
to the (approximate) conservation of the adiabatic charge, we follow the evolution of the 
/-balls while the coefficient of the logarithmic potential K is varied adiabatically. If the 
adiabatic invariance indeed guarantees the stability of the /-balls, the /-ball conhguration 
will be gradually deformed into a Gaussian prohle with a different value of K, while the 
adiabatic charge I is conserved. 

We introduce the time variation of K as 


Kit) 


Kq 

(1 -|- amt) ’ 


( 64 ) 


where Kq is the initial value of // at t = 0, and a is the coefficient of the time variation. 
For a 1, K varies much more slowly than the oscillation period, and therefore, the 
/-ball is expected to evolve into a Gaussian prohle with a different value of K. Thus, we 
expect that the /-ball radius R and R 1/2 evolve as 


\ Kq m 


/ln2 , /9 1 

Ri/2 = \ ^(1 + amtfl^ — . 

\ Kq m 

The typical time scale over which the radius R changes signihcantly is 


I 


-1 


am 


2 1 -|- amt 


-1 


am 


( 65 ) 

( 66 ) 


( 67 ) 


Therefore we need to follow the evolution of the /-balls for a sufficiently long period 
(3> Atfi) in order to see the adiabatic deformation. 

How small should a be for the /-ball deformation to be adiabatic? To answer this 
question, let us consider the deformation induced by excitations of the wave packets 
inside the /-ball. The typical time scale for the wave packet to transverse the entire 
region of the /-ball can be estimated as 


At50 


R 2 
_ _ 

Vg ~ KQm' 


( 68 ) 
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where Vg is the group velocity Vg = du/dk{k ~ 1/-R) — Kq/2. For adiabatic defor¬ 

mation of the /-ball, this propagation scale should be much smaller than At^, i.e., 
At 5 ^ AtR, which constrains a and Kq as 

a Kq. (69) 


If this condition (69) is met, the /-ball would deform adiabatically. 

The adiabatic charge I of the /-ball is expected to be conserved during the adiabatic 
deformation. 


RoKo 






(70) 


where the subscript 0 means that the variable is evaluated at f = 0 (see Eq. (54)). As 
long as // -C 1, the oscillation frequency is given by m up to a correction of order K, and 


so. 




+ 0{K). 


Therefore, the oscillation amplitude at the center, <I)c, should evolve with time as 

-DI2 


ic{t) ~ 4c,o (yyyj = 


( 71 ) 


(72) 


up to a small correction of order K. With this approximation, the effective amplitude <I>c 
evolves similarly, $c(^) ~ ‘^*c(^) (see Eq. (61)). 

First let us show the results for the case of D = 1, where we set the box size L and 
the number of grid N to 


L = 450/m(t = 0), N = 2048, At = lO"^ —. (73) 

m 

We have followed the evolution of the /-ball from f = 0 to f = 10^/m for Kq = 10“^ and 
a = 10“^. As a result the coefficient K{t) evolves from Kq to (approximately) iFo/lOO, 
and the /-ball radius is expected to become larger by a factor of 10. In Fig. 2, we show the 
numerical results. The top two panels show snapshots of the spatial distribution of p at 
t = 0, 10^, 10^, 10^ [1/m] with linear and logarithmic scales. One can see that the radius 
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FIG. 2: Numerical results for a = 10“^ and Kq = 10“^ in the case of D = 1. The top panels 
show the snapshots of the spatial distribution of p at t = 0,10^, 10^, 10^[l/m] in linear and 
logarithmic scales. The bottom panels show the evolution of R 1/2 and <hc from t = 10/m to 
10^[l/m]. The red (green) line shows the numerical (analytical) result (See Eqs. (66) and (72)). 


of J-ball becomes larger and its amplitude at the center becomes smaller as expected. 
The two bottom panels show the time evolution of R 1/2 and <hc in a very good agreement 
with the analytic estimate. This result clearly shows that the adiabatic charge I of the 
J-ball is indeed conserved, and that the J-ball conhguration follows the analytic solution 
obtained as the minimal energy state for a given adiabatic charge. 

We have similarly studied the deformation of the J-balls in the case of iJ = 2 and 
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D = 3 for i^o = 10 ^ and a = 10 We set the grid number and the box size as 


(74) 

(75) 


N = 256^ L = 100/m, for D = 2 
N = 128^ L = 50/m, for D = 3 

and At = 10~^/m for both cases, and followed the evolution from f = 0 to f = 10^/m. 
The results of the simulations are summarized in Fig. 3. From the top panels, we can 
see that as the coefficient K becomes smaller, the /-ball radius becomes larger. This 
deformation follows the analytic solutions obtained under the assumption of the adiabatic 
charge conservation, as can be see from the middle and bottom panels in Fig. 3. 

We have conhrmed the adiabatic deformation of /-balls for a = 10“^. For a larger 
a, however, the deformation of /-ball is no longer adiabatic (see (69)), and it does not 
follow the analytic prohle as the adiabatic charge is not conserved. In Fig. 4, we show 
the results of the case of D = 1 with {Kq, a) = (10“^, 10“^), for which the condition (69) 
is (marginally) broken. The /-ball does not have much time to deform itself in response 
to the change of K. As one can see from Fig. 4, the /-ball conhguration does not follow 
the Gaussian prohle any more and the evolution of the radius and the amplitude do not 
match the analytic one. 

V. DISCUSSION AND CONCLUSIONS 

The longevity of oscillons//-balls was a puzzle, and it was conjectured in Ref. [3] that 
it is due to conservation of the adiabatic charge /, in much the same way as the Q-balls 
are stable due to conservation of the global U(l) charge Q. There, it was numerically 
conhrmed that the /-ball held conhguration for the LR mass term potential agrees very 
well with the analytic one derived as the lowest energy state for a given adiabatic charge, 
giving a strong support for the conjecture. 

In this paper we have hrst proved the adiabatic current conservation for a potential that 
allows periodic motion with a separable form (11) even in the presence of non-negligible 
gradient energy. We have found that the modihed energy density p is a constant of motion, 
which plays a crucial role in the proof. Then we have determined a possible form of the 
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scalar potential that allows periodic motion to be the LR mass term potential (40). We 
have derived the J-ball solution using the Gaussian ansatz [3], under the condition that 
the adiabatic charge conservation. Finally we have numerically conhrmed the evolution 
and properties of the /-balls. In particular, we have followed the adiabatic deformation 
of the /-balls while the coefficient of the logarithmic potential K varies sufficiently slowly 
with time, and con&med that the numerical results perfectly agree with the analytic 
estimates based on the adiabatic charge conservation. We have also checked that, once 
the adiabatic condition (69) is violated, the /-ball starts to spread out and its evolution 
does not follow the analytic estimate.® Thus, our results unambiguously show that the 
stability (longevity) of the oscillons//-balls is due to the (approximate) conservation of 
the adiabatic charge. 

We have followed the evolution of the /-ball starting from the Gaussian solution. As 
already shown in Ref. [3], /-balls with the Gaussian prohle are formed if one started with 
the spatial homogeneous initial condition. For K > 0, there is an instability band, whose 
growth rate is of order mK. Therefore the adiabatic condition is only marginally satished. 
It is of interest to study the adiabatic charge conservation during the /-ball formation, 
and we leave it for future work. 

We have focused on the LR mass term potential which enables the periodic motion, 
making the /-ball absolutely stable at the classical level. For the other types of potential 
dominated by the quadratic potential, the adiabatic charge is approximately conserved. 
The oscillons//-balls for such potentials are considered to be long-lived due to the ap¬ 
proximate conservation of the adiabatic charge. Let us denote the deviation from the 
LR mass term potential by a small parameter e. The scalar dynamics is no longer given 
by the separate form (34) because of the deviation. In particular the trajectory over one 
period is not closed by an amount of e. Noting that the adiabatic invariant in the classical 
mechanics is a well conserved quantity, and its variation is exponentially suppressed for a 
small breaking of the adiabaticity [22], it is plausible that the approximate conservation of 

® If iy is increased rapidly with time, on the other hand, the /-balls are considered to split into smaller 
pieces. 
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the adiabatic charge accounts for the longevity of oscillons observed in various numerical 
simulations [8-10]. Intriguingly, it was shown that the oscillons in the small-amplitude 
regime emit scalar waves in an exponentially suppressed way [4], The violation of the 
adiabatic charge may enable us to understand the lifetime of the /-ball analytically. 

The LR mass term potential (and other types of potentials dominated by the mass 
term) appears in various cases. For instance, a scalar potential for flat directions in 
supersymmetric theory is often approximated by such LR mass term potential. The /-balls 
may be formed in the early Universe, and they may play an important cosmological role, 
especially if they are sufficiently long-lived. In the end of the day the /-ball may decay 
by violation of the adiabatic charge or quantum processes. We leave the cosmological 
application of the /-balls for future work. 
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FIG. 3: Numerical results for a = 10“^ and Kq = 10“^ in the cases oi D = 2 and D = 3. 
The left (right) panels show the results for ZD = 2 (3) case. The top two panels show snapshots 
of the spatial distribution of p at t = 0,10^, 10^[l/m], where Nx,Ny and represent the grid 
point number of the lattice. The middle (bottom) panels show the evolution of R 1/2 (d>c) from 
t = 1/m to 10^/m in a very good agreement with the analytical ones (66) and (72). 
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FIG. 4: Same as Fig. 2 but for a = 10 The /-ball deformation is no longer adiabatic, and 
it does not follow the analytic one. 
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